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ABSTRACT 

The CM Draconis system comprises two eclipsing mid-M dwarfs of nearly equal mass in a 1.27-day 
orbit. This well-studied eclipsing binary has often been used for benchmark tests of stellar models, 
since its components are amongst the lowest mass stars with well-measured masses and radii (< 1% 
relative precision). However, as with many other low-mass stars, non- magnetic models have been 
unable to match the observed radii and effective temperatures for CM Dra at the 5-10% level. To 
date, the uncertain metallicity of the system has complicated comparison of theoretical isochroncs 
with observations. In this Letter, we use data from the SpcX instrument on the NASA Infrared 
Telescope Facility (IRTF) to measure the metallicity of the system during primary and secondary 
eclipses, as well as out of eclipse, based on an empirical metallicity calibration in the H and K near- 
infrared (NIR) bands. We derive a [Fe/H] = —0.30 ± 0.12 that is consistent across all orbital phases. 
The determination of [Fe/H] for this system constrains a key dimension of parameter space when 
attempting to reconcile model isochrone predictions and observations. 



1. INTRODUCTION 

CM Draconis (GJ 630.1 AC, Gliese & Jahreiss 1991, 
hereafter CM Dra) is an eclipsing binary system that con- 
sists of a pair of very similar M dwarfs. The eclipses are 
nearly total (89% and 99.9% of light is blocked from the 
primary and secondary, respectively). The first spectro- 
scopic mass measurements were made by Lacy (1977), 
who noted that the two stars did not lie near the the- 
oretical zero-age main sequence for Population I stars. 
A more recent analysis by Morales et al. (2009) has de- 
termined a primary mass Mi = 0.2310 ± 0.0009 Mq and 
radius i?i = 0.2534 ± 0.0019 i?0, with a secondary mass 
M2 = 0.2141 ± 0.0010 Mq and radius R2 = 0.2396 ± 
0.0015 i?0. The system's appreciable proper motion, rel- 
atively large radial velocity, and the existence of a com- 
mon proper motion white dwarf companion have been 
claimed as evidence that the system may be a Popula- 
tion II member (Lacy 1977). 

Observationally determined radii and effective temper- 
atures of low- mass stars (M < 0.8 Mq) disagree signifi- 
cantly with theoretical isochrones (Lopez-Morales 2007; 
Morales et al. 2008; Torres et al. 2010; Stassun et al. 
2012), and CM Dra is no exception (Morales et al. 2009; 
Feiden et al. 2011). This disagreement was generally 
at the 5% - 15% level, but a more recent study by 
(Feiden & Chaboyer 2012) uses a finely sampled grid 
of stellar parameters and an improved equation of state 
compared to previous studies, reducing the disagreement 
to the sub-5% level for a majority of systems. Magnetic 
inhibition of convective efficiency or an increase of cool 
starspots at polar latitudes are thought to be the pri- 
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mary cause (Chabricr et al. 2007; Morales et al. 2010; 
Feiden & Chaboyer 2012) of this discrepancy, but un- 
known parameters such as age and metallicity make it 
difficult to compare observations with theory. In the case 
of CM Dra, models attempting to mimic the interactions 
of a magnetic field with the stellar interior can be made 
to match the observations, but often require a super-solar 
metallicity, a large starspot coverage, or very strong mag- 
netic fields (Morales et al. 2010; MacDonald & MuUan 
2012). 

The metallicity of the CM Dra system merits special 
consideration because it can constrain one dimension of 
the modeling problem for a system in the fuUy-convectivc 
regime (< 0.35 Mq) where there are only four other well- 
characterized systems: KOI-126 (Carter et al. 2011), 
Kepler 16 (Doyle et al. 2011; Bender et al. 2012), Ke- 
pler 38 (Orosz et al. 2012), and LSPM J1112+7626 
(Irwin et al. 2011). In the first three cases, metallicity 
estimates can be obtained from their higher-mass pri- 
mary stars, but no such primary is available for CM 
Dra, and so metallicities must be determined from the 
low-mass stars alone. 

Attempts to measure the metallicity of CM Dra using 
synthetic spectral templates in both the optical and near- 
infrared (NIR) have resulted in metallicity values ranging 
from solar (Gizis 1997) to [M/H] = -1.0 (Viti et al. 
2002). However, such techniques rely on model atmo- 
spheres of cool stars that have incomplete line lists and 
significant uncertainties in the complicated molecular 
chemistry occurring in the outer atmospheres. An alter- 
native approach is to use an empirical relationship based 
on equivalent widths of specific, narrow regions of the 
spectrum (Rojas-Ayala et al. 2010, 2012; Terrien et al. 
2012). 

In this Letter, we utilize such an empirical relation to 
derive the [Fc/H] of CM Dra from spectra taken at a vari- 
ety of orbital phases: primary eclipse, secondary eclipse, 
and out-of-eclipse. By using a technique that does not 
rely heavily on low-mass stellar atmosphere models, and 
obtaining observations of the near-total eclipses that iso- 
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late the light from one star, we can obtain a reliable 
estimate of CM Dra's metallicity. In Section 2 we de- 
scribe our SpeX data and the empirical relation used to 
derive metallicitics. In Section 3 we derive our metallic- 
ity estimate and compare it to values in the literature. In 
Section 4 we discuss the impact of this metallicity con- 
straint on stellar models of the CM Dra system. Finally, 
in Section 5 we summarize our results. 

2. OBSERVATIONS WITH THE IRTF-SPEX 
SPECTROGRAPH 

During May 2012, we observed CM Dra with the NASA 
Infrared Telescope Facility (IRTF) SpeX spectrograph 
(Rayner et al. 2003) as a part of a larger program to 
estimate metallicitics for several hundred M dwarfs. As 
part of this larger program, we have developed indepen- 
dent calibrations of metallicity ([Fe/H]) for the K and 
i?-bands of IRTF-SpeX M dwarf spectra (Terrien et al. 
2012). We operated SpeX in the short cross-dispersed 
(SXD) mode with the 0.3 x 5" slit, which produces R ~ 
2000 spectra from 0.8 - 2.4 ^m. We extract these spec- 
tra with the facility SpeXTool package (Cushing et al. 
2004). We telluric-correct and flux-calibratc the spectra 
using the xtellcor program (Vacca et al. 2003) and ob- 
servations of an AOV star (or similar) , at an airmass and 
observation time within 0.1 (sec z) and 1 hour of the tar- 
get, respectively. Our analysis of the reduced spectra is 
identical to that presented in Terrien et al. (2012), but 
here we note a few important points. 

Our K and i?-band metallicity relations are based on 
the equivalent widths of Na i and Ca i in the i^-band 
and K i and Ca i in the i7-band. The windows used 
to calculate the strengths of these features arc defined 
and explained in Terrien et al. (2012). These windows 
are shifted for each observation by a radial velocity (RV) 
offset derived using a K-hand spectral template (a high 
S/N spectrum of HD 36395 from Rayner et al. 2009). ^ 
The measured equivalent widths arc not sensitive to ro- 
tational broadening for a typical M dwarf rotational ve- 
locity (< 30 km/s), as the rotational broadening is small 
compared to the spectrograph resolution (R ~ 2000). For 
CM Dra in particular. Morales et al. (2009) used visible 
light echelle spectra with a spectroscopic resolution of 
^ 35000 to measure wsini sa 10 km/s for both com- 
ponents. We have simulated the impact of this v sini on 
the Terrien et al. (2012) windows, and find it affects our 
measurements of equivalent width and metallicity at less 
than 0.1%. 

The metallicity relations arc calibrated using a set of 
22 M dwarf companions to higher-mass stars that have 
well-determined metallicitics from the SPOCS catalog 
(Valenti & Fischer 2005), and yield [Fe/H], which can 
be used as a proxy for [M/H]. We have demonstrated that 
the two independent calibrations are consistent with one 
another (Terrien et al. 2012). Since the features in the 
H-hasaA are weaker than those used in the iiT-band, the 
iJ-band calibration requires S/N « 150 pixel" ^ while the 
K-hand requires S/N « 100pixer\ Our adopted error 
in the calibrations is ±0.12 dex, based on the the scatter 
of the metallicity estimates compared to the metallicity 
measurements of the primaries. 

^ http://irtfweb.ifa.hawaii.edu/~spcx/IRTF_Spcctral_Library/ 
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Figure 1. RV model (using RVLIN, Wright & Howard 2009) 
of CM Dra using orbital parameters from Morales ct al. (2009), 
showing the primary (solid line), secondary (dashed line), and 
gamma velocity (dotted line). The gray bars mark the phases 
at which our IRTF spectra were obtained: primary eclipse (Pri.), 
secondary eclipse (Sec), and out-of-eclipse (OoE). 

We observed CM Dra over a period of several days, dis- 
tributing observations over many phases, including one 
primary and two secondary eclipses, as well as two out-of- 
eclipse observations. Both of our out-of-eclipse observa- 
tions occurred close to maximum RV separation (Figure 
1), have high S/N (> 300), and contain full light from 
both components of CM Dra. We planned the three in- 
eclipse observations using the orbital parameters from 
Morales et al. (2009). We obtained coverage for 30-40 
minutes, centered on the times of maximum eclipse, by 
taking repeated 30 second exposures and combining them 
into ~ 2 minute exposures during extraction. Both pri- 
mary and secondary eclipses for CM Dra last a total of 
^ 1.5 hours. At the midpoint of the primary eclipse, the 
primary is ~ 89% obscured; at the midpoint of the sec- 
ondary eclipse, the secondary is ^ 99.9% obscured. We 
can empirically estimate the effect of spectral contam- 
ination from either component on our metallicity mea- 
surements by observing during both eclipses and out of 
eclipse. 

Around the time of maximum secondary eclipse, the 
secondary remains more than 98% obscured for approx- 
imately three minutes, and we obtained the S/N > 150 
required to perform both the K and iJ-band metallicity 
estimates. Since our calibration is valid only for single 
stars, the isolated light from the primary allows us to es- 
timate the metallicity of this component of the CM Dra 
system. Assuming the stars are coeval, this then provides 
our best metallicity estimate for the CM Dra system. 

3. RESULTS AND DISCUSSION 

Table 1 lists our individual CM Dra metallicity mea- 
surements and Table 2 summarizes them by epoch. The 
combined K and 7?-band measurements for the sec- 
ondary eclipse (effectively the metallicity estimate of the 
primary), as split into 2-minute exposures, have a mean 
of [Fe/H]= —0.30 dex with a spread of cr = 0.03 dex, 
which is much smaller than the absolute uncertainty in 
our metallicity calibration relationship (Section 2). We 
therefore adopt the uncertainty in our calibration as the 
error for the metallicity estimate for CM Dra. We find 
CM Dra to be slightly metal-poor relative to the Sun, 
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with the caveat that we are not sensitive to abundance 
patterns which could yield differences between [Fe/H] 
and [M/H]. 

Interestingly, the primary eclipse and out-of-eclipse 
metallicity estimates do not deviate significantly from the 
measurement of the secondary eclipse. Since these stars 
are closely bound and are of similar mass (7% difference) 
and spectral type (both dM4.5), this supports the valid- 
ity of the assumption that the two stars are coeval, since 
differing compositions would likely yield different feature 
strengths between individual and combined light spectra. 

The differences that do exist between the secondary 
eclipse observations and the primary /out-of-eclipse ob- 
servations are only present in the iJ-band estimate, and 
are strongly correlated with airmass: when airmass is 
used to linearly model [Fe/Hjn-band, the multiple cor- 
relation coefficient ~ 0.78. These differences are 
therefore likely byproducts of the sensitivity of the H- 
band technique to atmospheric contamination, as noted 
in Terrien et al. (2012). In addition to the sensitivity of 
the K I feature noted in the previous paper, in perform- 
ing the CM Dra analysis we have also discovered a signif- 
icant sensitivity of the Ca i feature to a step in the xtell- 
cor routine. Specifically, the user must manually scale 
the strengths of Brackett lines in a model spectrum of 
Vega in order to remove them from the observed telluric 
standard. It is not immediately obvious how to correctly 
scale a particular Brackett line in the vicinity of both the 
iJ-band Ca i feature and a telluric absorption feature. 
Repeated iterations of our analysis allowed us to develop 
a consistent scaling based on requiring a certain smooth- 
ness at the edge of a telluric CO absorption feature across 
the multiple telluric standard observations for CM Dra. 
We are studying methods of dealing with the dependency 
of the H-hand technique on telluric contamination us- 
ing high-resolution if-band spectra from the SDSS-III 
APOGEE instrument (AUende Prieto et al. 2008). 

Since CM Dra has both a high absolute RV 
(« 100 km s~^) and high orbital RV amplitude (w 
70 km s~^), it is important to ensure that the windows 
used to calculate feature strengths are shifted appropri- 
ately in RV (the window widths remain the same). Fig- 
ure 2 shows the May 7 out-of-eclipse observations (near 
maximum RV separation), the feature windows shifted 
by our analysis routine, and an RV-shifted BT-Settl spec- 
trum for reference (AUard et al. 2003, 2007, 2011). 

Our estimate of the metallicity of CM Dra is the first 
in a long series of composition studies to converge on a 
precise answer. Comparison of individual features with 
models in both the optical and NIR (e.g. atomic fea- 
tures and a CO band, Viti et al. 1997, 2002) have con- 
cluded that CM Dra is metal poor ([M/H] « -1.0), which 
would be consistent with it being a halo/population II 
object. However, similar studies of other features (e.g. 
CaH and TiO, Gizis 1997) and IR colors (Leggctt 1998) 
have preferred a near-solar composition ([M/H] « 0.0), 
and multiple studies (Viti et al. 1997; Leggett 1998) 
have noted difficulty in reconciling metallicity determi- 
nations for different wavelength regions. Recent studies 
have even included attempts at detailed line modeling 
(Kuznetsov et al. 2012). Our finding that CM Dra is 
slightly metal poor ([Fc/H] « —0.30) lies between the 
aforementioned studies that find a range of mctallici- 
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Figure 2. Shifted spectral windows (from Terrien et al. (2012)) 
used to calculate the EWs (grey regions) compared to a pair of 
BT Settl models of T^a = 3100 K, logg = 5.0 and [M/H] = -0.5, 
Doppler-shifted to the RVs expected at the time of the first out- 
of-eclipse observation and convolved to R = 60,000 (black). The 
red, dashed line is the observed R = 2000 IRTF spectrum for 
comparison. 



ties between metal poor and solar, and is not consistent 
with studies that find the system to be metal poor. Un- 
like these previous studies, our empirical estimate of the 
metallicity of CM Dra does not suffer from the difficulties 
of applying low-mass stellar atmospheric models, such as 
missing or poorly defined opacities. There is one pub- 
lished metallicity estimate using a similar technique to 
ours, from Rojas-Ayala et al. (2012), which finds [Fc/H] 
= —0.39 ± 0.17 dcx using features in the X-band. This 
estimate is consistent with ours, but to the best of our 
knowledge is not timed for a specific epoch and so likely 
contains mixed light. 

4. IMPACT ON STELLAR MODEL INTERPRETATION 

Since CM Dra is among the lowest mass systems with 
well-determined stellar parameters, it is an important 
benchmark for low-mass stellar models, and will play 
an important role in work to reconcile the discrepancy 
between low-mass models and observations. The most 
recent study to address stellar models and CM Dra, 
Feiden & Chaboyer (2012), still finds a discrepancy of 
^ 3% between the observed and modeled radius for both 
components. This was derived using a metallicity prior 
of —1 < [Fc/H] < 0, which allowed a best-fit metal- 
licity of as input to the Dartmouth Stellar Evolution 
Program (DSEP, Dottcr et al. 2007, 2008; Feiden et al. 
2011; Feiden & Chaboyer 2012). Figure 3 A shows a 4 
Gyr DSEP isochronc with [Fe/H] = -0.3. The model radii 
are i?A, model = 0.2374 Rq and i?B,modci = 0.2213 Rq. 
Using the radii measured by Morales et al. (2009), the 
relative errors, defined as 



SR 

Robs 



Robs — R 



model 



Rob 



(1) 



are SRA/Robs = 0.063 and SRs/Robs = 0.076. If 
we take instead the l-cr upper limit for metallicity of 
[Fc/H]= —0.18, we find that the discrepancies decrease 
to 5RA/Robs = 0.052 and SRs/Robs = 0.065. Higher 
metallicitics yield larger radii (somewhat ameliorating 
the problem, as in e.g. Spada & Demarquc 2012), so the 
true extent of the discrepancy for CM Dra was almost 
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Table 1 

CM Dra Metallicity Measurements 



/^-band //-band 



RJD^ 


In Eclipse? 


Phase'' 


EWn;, \A] 

' ' • • IN a *-j 


EWr„ [Al 


S/Nk 


fFe/HlK 


EWpa \K\ 


eWk f Al 


S/Nh 


fFe/HlH 


56055.03708 


No 


0.229 


3.58 ± 0.07 


2.28 ± 0.12 


312 


-0.30 


1.25 ± 0.11 


0.57 ± 0.05 


384 


-0.40 


56056.90978 


No 


0.706 


3.44 ± 0.06 


2.20 ± 0.07 


344 


-0.33 


1.38 ± 0.04 


0.63 ± 0.02 


404 


-0.37 


56050.93011 


Primary 


0.992 


3.45 ± 0.08 


2.27 ± 0.10 


246 


-0.32 


1.06 ± 0.07 


0.78 ± 0.04 


274 


-0.38 


56050.93380 


Primary 


0.994 


3.46 ± 0.07 


2.25 ± 0.09 


260 


-0.32 


1.06 ± 0.07 


0.81 ± 0.03 


285 


-0.36 


56050.93749 


Primary 


0.997 


3.42 ± 0.08 


2.23 ± 0.10 


243 


-0.33 


1.12 ± 0.06 


0.84 ± 0.03 


264 


-0.34 


56050.94118 


Primary 


0.000 


3.40 ± 0.07 


1.97 ± 0.10 


235 


-0.35 


0.95 ± 0.06 


0.77 ± 0.04 


257 


-0.38 


56050.94486 


Primary 


0.003 


3.16 ± 0.08 


2.23 ± 0.08 


245 


-0.36 


0.97 ± 0.08 


0.79 ± 0.03 


268 


-0.36 


56050.94855 


Primary 


0.006 


3.39 ± 0.07 


2.31 ± 0.08 


277 


-0.32 


1.03 ± 0.05 


0.78 ± 0.03 


305 


-0.37 


56050.95316 


Primary 


0.010 


3.53 ± 0.06 


2.21 ± 0.06 


368 


-0.31 


1.03 ± 0.04 


0.81 ± 0.02 


406 


-0.35 


56052.83425 


Secondary 


0.493 


3.49 ± 0.10 


2.20 ± 0.13 


162 


-0.32 


1.27 ± 0.09 


0.84 ± 0.04 


189 


-0.28 


56052.83794 


Secondary 


0.496 


3.36 ± 0.11 


2.26 ± 0.12 


171 


-0.33 


1.33 ± 0.08 


0.84 ± 0.04 


195 


-0.23 


56052.84162 


Secondary 


0.499 


3.81 ± 0.10 


2.39 ± 0.12 


168 


-0.26 


1.41 ± 0.08 


0.84 ± 0.04 


193 


-0.25 


56052.84531 


Secondary 


0.501 


3.53 ± 0.12 


2.34 ± 0.11 


184 


-0.30 


1.23 ± 0.07 


0.85 ± 0.03 


210 


-0.27 


56052.84992 


Secondary 


0.505 


3.34 ± 0.09 


2.22 ± 0.09 


231 


-0.34 


1.17 ± 0.06 


0.90 ± 0.03 


261 


-0.26 


56052.86108 


Secondary 


0.514 


3.40 ± 0.09 


2.40 ± 0.12 


194 


-0.32 


0.96 ± 0.06 


0.85 ± 0.03 


225 


-0.32 


56052.86477 


Secondary 


0.517 


3.27 ± 0.08 


2.29 ± 0.10 


238 


-0.34 


0.98 ± 0.05 


0.87 ± 0.03 


268 


-0.30 


56054.10690 


Secondary 


0.496 


3.51 ± 0.10 


2.24 ± 0.12 


172 


-0.31 


1.04 ± 0.08 


0.89 ± 0.04 


197 


-0.30 


56054.10920 


Secondary 


0.498 


3.41 ± 0.10 


2.48 ± 0.11 


177 


-0.31 


1.19 ± 0.08 


0.81 ± 0.04 


202 


-0.33 


56054.11150 


Secondary 


0.500 


3.55 ± 0.11 


2.11 ± 0.13 


172 


-0.32 


1.21 ± 0.08 


0.83 ± 0.05 


197 


-0.27 


56054.11381 


Secondary 


0.502 


3.41 ± 0.10 


2.35 ± 0.12 


186 


-0.32 


1.18 ± 0.07 


0.87 ± 0.03 


213 


-0.31 


56054.11611 


Secondary 


0.503 


3.42 ± 0.33 


2.37 ± 0.24 


193 


-0.31 


1.18 ± 0.11 


0.82 ± 0.17 


217 


-0.32 


56054.11842 


Secondary 


0.505 


3.46 ± 0.10 


2.31 ± 0.12 


187 


-0.31 


1.18 ± 0.07 


0.84 ± 0.04 


215 


-0.32 


56054.12129 


Secondary 


0.507 


3.35 ± 0.08 


2.28 ± 0.11 


211 


-0.33 


1.10 ± 0.07 


0.91 ± 0.03 


240 


-0.31 


56054.12389 


Secondary 


0.509 


3.55 ± 0.09 


2.52 ± 0.10 


239 


-0.28 


1.27 ± 0.06 


0.86 ± 0.05 


271 


-0.28 



Average RJD of combined exposure. 
^ Phase calculated using elements from Morales et al. (2009). 



Table 2 

Averaged CM Dra Metallicity Measurements per Epoch 



Epoch 


RJD 


[Fc/HJk 




[Fc/HIh 




Out of eclipse 1 


56055 


-0.30 




-0.40 




Out of eclipse 2 


56056 


-0.33 




-0.37 




Primary Eclipse 


56050 


-0.33 


0.02 


-0.36 


0.02 


Secondary Eclipse 1 


56052 


-0.32 


0.03 


-0.27 


0.03 


Secondary Eclipse 2 


56054 


-0.31 


0.01 


-0.31 


0.02 



Note. — The summary by epoch of the metallicity mea- 
surements performed on CM Dra. The metallicity measure- 
ments are the means and standard deviations calculated for 
each epoch of observation. 

certainly underestimated in Feiden & Chaboyer (2012) 
due to the relatively high allowed metallicity of 0. 

Besides metallicity, another significant constraint on 
the stellar models is age. The age of the CM Dra sys- 
tem is supposedly well-constrained by the cooling age 
of its white dwarf companion (Bergeron et al. 2001; 
Morales et al. 2009). Figure 3 B shows the evolution- 
ary tracks for models with the masses of CM Dra A/B 
and reveals the extent to which age can affect the radius 
discrepancy for this system. With [Fe/H] = — 0.30±0.12, 
CM Dra B can only be realistically fit with ages of 100- 
200 Myr. However, the estimated 4 Gyr age of the white 
dwarf companion and the high space velocity of the sys- 
tem suggest that this system is much older. We conclude 
that with constraints on both age and metallicity, the 
non-magnetic models for this system remain discrepant 
at the 5-7% level. 

With the present constraint on metallicity and a 
long history of photometric and spectroscopic moni- 
toring and modeling, CM Dra provides a strong test 



of low-mass stellar models. KOT126 (Carter ct al. 
2011; Feiden et al. 2011), Kepler 16 (Doyle ct al. 2011; 
Winn ct al. 2011) and LSPM J1112+7626 (Irwin ct al. 
2011) provide constraints in the fuUy-convective region 
as well, as will the recently-discovered Kepler 38 B 
(Orosz et al. 2012), when a larger number of observa- 
tions are able to constrain its mass and radius to the 
~ 1% level. The discovery papers for both KOT126 and 
Kepler 16 (Carter et al. 2011; Doyle et al. 2011) pro- 
vide spectroscopic estimates of metallicity of the higher- 
mass components of the systems, and for KOT126 these 
are supplemented by a model-derived age constraint from 
Feiden et al. (2011). Ahhough KOT126 is well-fit by the 
appropriate stellar models, CM Dra and Kepler 16 still 
deviate significantly (5-7% and 3%, respectively). Kepler 
16 is of particular note due to its well-constrained mass 
which has been independently spectroscopically verified 
(Bender et al. 2012). LSPM J1112-K7626B still devi- 
ates at the 3% level as well, although its metallicity 
and age are not as well constrained (Irwin et al. 2011; 
Feiden & Chaboyer 2012). As observations of these ob- 
jects tighten the constraints on stellar models, the dis- 
crepancy between models and observations persists, pos- 
sibly pointing the way toward new physics that must be 
included in the models, as mentioned in Section 1. Im- 
portantly, CM Dra is one of only two known systems 
with well-constrained masses and radii (along with KOT 
126) that are composed of low-mass stars in very short 
(« 1 day) periods, and that will be able to provide di- 
rect tests of models which explain the observed inflated 
radii of low-mass stars with magnetic fields produced by 
rotationally-powercd dynamos. A thorough understand- 
ing of low-mass stars, and especially well-constrained 
mass-radius relations, will also be particularly important 
for M dwarf planet searches (Nutzman & Charbonncau 
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Figure 3. Stellar models from the Dartmouth Stellar Evolution 
Program (Feiden & Chaboyer 2012). (A) Mass-radius diagram 
of stars that have 0.2 < M < 0.25 Mq and precise masses and 
radii. CM Dra A/B are shown in black, along with the DSEP 
isochrones for an age of 4 Gyr and our estimated metallicity val- 
ues [Fe/H]= -0.30 ± 0.12. Also plotted in gray are Kepler 16B 
and KOI-126B/C, along with their respective best-fit isochrones, 
from Feiden & Chaboyer (2012). The recently-discovered Kepler 
38 B (Orosz et al. 2012), which will be a strong test of isochrones 
when its stellar parameters are more tightly constrained, is also 
shown. (B) The evolutionary tracks for both components of CM 
Dra with the measured metallicity {[Fe/H] = —0.3), along with the 
same tracks for our l-cr upper limit metallicity estimate ([Fe/H] 
= -0.18). The tracks for the l-cr lower limit ([Fe/H] = -0.42) are 
indistinguishable from those for [Fe/H] = —0.30 in this plot. 

2008; JVIahadevan et al. 2010; Quirrenbach et al. 2010), 
as radius is a crucial input in the calculation of the loca- 
tion of the habitable zone (e.g. Boyajian et al. 2012). 

5. CONCLUSION 

We have estimated the metallicity of the CM Dra sys- 
tem using two empirical NIR relations applied to observa- 
tions in and out of eclipse. We find it to be slightly metal 
poor, with [Fc/H]= — 0.30±0.12 dex. This measurement 
is consistent with estimates obtained out of eclipse and 
during primary eclipse, suggesting that the components 
are very similar in composition, since they are close in 
mass and spectral type. When used as input to DSEP 
stellar models, this study and studies of the similar sys- 
tems KOI-126, Kepler 16, and LSPIVI J1112+7626 show 
that the discrepancies between measured and predicted 
radiii persist at the 3-7% level. The precise determina- 
tion of its metallicity, as well as the constraint on its age 
from the proposed white dwarf companion, enables CIVI 
Dra to be one of the best test cases for low-mass stel- 
lar models. Further improvement in the understanding 



of low-mass stars will require both observational work in 
finding and characterizing low-mass stars, and theoreti- 
cal work in modeling the additional physics necessary to 
match the observed cases. 
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